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1. Introduction 

In recent years, many efforts have been devoted to the definition of discrete analogue of con¬ 
tinuous methods and objects, like variational calculus, Lagrangian and Hamiltonian systems, 
etc. We refer to [9] and m for an overview of this subject. Although different points of 
view can be adopted, a common point of all these works is to understand in which extend a 
discrete object will differ from the original continuous one, as for example symmetries, first 
integrals and related properties. 

Following our previous work [5J, we continue to investigate the discrete calculus of variations 
in the formulation given for example in [5] following the classical work of m and [91. This 
formulation, which is obtained using the formalism of discrete embedding (see j6]), has the 
property to mimic very closely the continuous case. As an example of this difference in the 
formulation, on a can compare the discrete Euler-Lagrange equation given in [9] and the one 
given in [5] or [6) . 

In [3], we have derived the discrete Helmholtz’s condition for second order difference equa¬ 
tions, i.e. conditions under which such equations can be written as a discrete Euler-Lagrange 
equation. In this paper, we want to discuss under which condition a set of difference equa¬ 
tions can be written as a discrete Hamiltonian system. This subject is studied in PQ and [7] 
where they obtain some conditions. However, the underlying framework concerning the varia¬ 
tional integrators is not easily related the continuous one, as well as the proofs which are given. 

In this paper, we first derive in the framework of discrete embedding, a coherent definition 
of discrete Hamiltonian systems. Comparison of our presentation with previous attempts, 
in particular ([HUM) will be given. Second, we derive the Helmholtz conditions in the 
Hamiltonian case following the usual self adjoin characterization of the differential operator 
associated to the Hamiltonian system as exposed for example by Santilli m 
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PART I 

NOTATIONS 


We remind several notations and results developed in |6j in the framework of the discrete 
embedding. We refer to this article for more details and proofs. 

Let N G N and a, b G R with a < b and let h = (b — a)/N. We denote by T the 
subspace of K defined by T = KL D [ a,b\ where /iZ = {hz\z G M}. The elements of T are 
denotes by t^ for k = 0, N. We denote by T + = T\{f/v}, T _ = T\{to} and = T + (~lT _ . 

We denote by C([a, 6],M d ) the set of functions x : [a, b\ R d , d G N* and by C*([a, b], M d ) 
the set of i-th differentiable functions. 

We denote C(T, M d ) the set of functions with value in over T and 6 q(T. the subset 
defined by 

CoiT^) = {G G C(T,R d ), G 0 = G n = 0}. (1) 

The discrete derivatives A and V defined over are respectively defined by the forward and 
backward finite differences operator with value in C'(T + ,M rf ) and C'(T - ,M rf ) respectively. 
The discrete antiderivative is denoted by Ja- 

We define the set to be the discrete functions F € C(T, U. d ) such that Ja(F * F) is 

well defined where a denotes the discrete product over C(T,M d ) defined in ([6],§. 3.3.1). 

We have the following results (see |6j, §.5.3.2 and §.6.2) : 

Theorem 1 (Discrete integration by part formula). - Let F € C(T,R d ) and G <E 
Cc^TT,JR 0 *), we have 

[J A (F*A{G))] N = — [Ja(F -k (VG))] n , 

[MF*V(G))] N = ~[MF*(AG))] N . [Z) 

Lemma 2. — Let F G C(T, M d ) such [Ja{F * G)] N = 0 for all G G C'o(T. M rf ) then F^ = 0 
for k = 1, N — 1. 


PART II 

DISCRETE HAMILTONIAN SYSTEMS 


Discrete Hamiltonian systems have already been defined by many authors as for example in 
(naisuni)- However, these formulations do not follow the strategy of discrete embedding 
which is aimed as providing an explicit relation between algebraic and analytic continuous 
structures and there discrete analogue (see As a consequence, we provide a self-contain 
derivation of discrete Hamiltonian systems in this setting. A comparison with previous ap¬ 
proaches is given. 
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2. Reminder about Hamiltonian systems 

2.1. Hamiltonian systems. — For simplicity we consider time independant Hamiltonian 
and Lagrangian. 


Definition 3 (Classical Hamiltonian). A classical Hamiltonian is a function H : R rf x 
such that for ( q,p ) £ C 1 ([a, b\, M d ) x C 1 ([a, b], M. d ) we have the time evolution of (q,p) 
given by the classical Hamilton’s equations 


q 

P 


_ 9H{q,p) 
dp 

_ dH(q,p) 
dq 


(3) 


T Y7H — (9H_ 6H _\T w j iere j denotes 


A vectorial notation is obtained posing 2 = ( q,py and VH = ( 75 ^- 575 ^ 
the transposition. The Hamilton’s equations are then written as 

dz 


dt 


= J ■ VH, 


(4) 


where J = 


0 Id 

—Id 0 


with Id the identity matrix on M. d denotes the symplectic matrix. 


We also denote by Xu the associated vector field defined by 

_dH dH 

Xu — o dq dp. 

oq op 


(5) 


An important property of Hamiltonian systems is that there solutions correspond to critical 
points of a given functional, i.e. follow from a variational principle. 


Theorem 4- The points ( q,p ) £ C 1 ([a, b], M d ) x C 1 ([a, b], M 0 ') satisfying Hamilton’s equa¬ 
tions are critical points of the functional 


C H : C , 1 ([o, 6 ],R d ) x C\[a,b\, 


( 6 ) 


where Lu ■ 


( q,P ) 1 —> T H {q,p)= f L H (q{t),p(t),q(t),p(t)) 

J a 

—> M is the Lagrangian defined by 


L H (x,y,v,w) = y-v-H(x,y). 


(7) 


2.2. Lagrangian versus Hamiltonian. — As we have already a consistent discrete 
theory of Lagrangian system, we will used this derivation of Hamiltonian systems in order to 
define discrete Hamiltonian systems. Doing so, we will see that all the objects are related by 
the discrete embedding procedure. 


Let L be a Lagrangian and denote by EL the corresponding Euler-Lagrange equation given 
by 
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dL 

We assume that the Lagrangian is admissible, i.e. that the map v —>• —— is invertible for all 

dv 

( q , t ) € M n x K. As a consequence, we can introduce the moment variable 

dL 

p=- 7 ^(q,q,t), (9) 

in order to rewrite the Euler-Lagrange equation as a first order system of differential equations 
given by 

q = g{q,p,t), 


dL 

p = -g^{q,g{q,p,t),t), 


( 10 ) 


dL 


where g is the inverse of ——. This change of variable will be called Legendre transform in the 

dv 

following. Introducing the Hamiltonian function 

H(q,p,t ) = L(q,g(q,p,t),t ) -pg(q,p,t), (11) 

one proves that equation (HOD is Hamiltonian with respect to H. 


3. Discrete Hamiltonian systems 

There exists many way to define discrete analogue of Hamiltonian systems. In this Section, 
we quickly remind some classical results which can be found for example in ram®) in a 
slightly different form. 

We follow the discrete embedding formalism exposed in |6,. Explicit comparisons with other 
formulations are given. In order to obtain a coherent definition of discrete Hamiltonian 
systems, we follow the usual continuous derivation, starting from a Lagrangian system to its 
Hamiltonian formulation via the Legendre transform. As we will see, all these structures can 
be obtained via a direct discrete embedding of the corresponding continuous structures. This 
property is exactly what is missed in the existing literature. 


3.1. Discrete Lagrangian systems. — The discrete analogue of the Euler-Lagrange equa¬ 
tion is given by (see [6]) : 


V 


dr 

—(Q,AQ,T) 


BL 

= —(Q,AQ,T). 


( 12 ) 


We introduce the set of variables 

dl 

P = —(Q,AQ,T), (13) 

which are the classical moment variables in classical mechanics. Thanks to the admissibility 
of L, we have A Q = g(P,Q,T). Using these variables, the discrete Euler-Lagrange equation 
is then equivalent to the discrete system 


P = -^-(Q,g(P,Q,T),T), 

dl 

VP = —(Q,g(P,Q,T),T). 


(14) 


The form of this discrete system is obtain from the continuous case taking care of the duality 
between the variable Q and P inducing a change between A and V for the discretisation of 
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the derivative. 


Introducing the Hamiltonian function 

H(Q, P, T) = L(Q, g(P, Q, T), T) - Pg(P Q, T ) (15) 

which is exactly the discrete form of the continuous Hamiltonian H associated to L. 

Remark 1. It must be noted that the previous result holds independently of the discrete 
embedding which is chosen. As a consequence, we have always the same Hamiltonian function. 
This is completely different from the result of Lall-West [5] where different Hamiltonian are 
introduced. 


Using the Hamiltonian, one obtain the Hamiltonian form of the Euler-Lagrange equation : 


A Q 
VP 


d H 

^8H 

dcy 


(16) 


3.2. Discrete Hamiltonian system. — The previous derivation of a discrete form of 
Hamiltonian systems using the discrete Euler-Lagrange equations leads to the following defi¬ 
nition : 


Definition 5 (Discrete Hamiltonian). - A discrete Hamiltonian is a function H : R d x 
R rf —> R such that for ( Q,P ) £ C^TT.R^) x C(T,R rf ) we have the time evolution of ( Q,P ) 
given by the discrete A—Hamilton equations (resp. V— Hamilton equations) 


(A Q) =§£, 

<yp) =-§> 

over . 


(17) 


Aa a convenient notation, in analogy with the continuous case, we introduce what we call a 
discrete vector field associated to the difference equation (1171) denoted by X a ,v,h and formally 
denoted by 

y dH i dH i nn 

X A.V,H = -Qpd A ,Q - -QQ d V,P- (18) 

With this notation, we have the following commutative diagram indicating the coherence of 
our construction with respect to the continuous one : 


The discrete diagram 


L h 


disc 


>H h 


dlap 


def 


^ Leqendre 

EL h -> X A y,H 
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where d.l.a.p denotes the discrete least action principle, is exactly the same as the continuous 
one 


L 


disc 




lap 


def 


EL 


Leqendre 

—-- > x H 


where l.a.p denotes the least action principle. All the objects are related via the discrete 
A -embedding as given by 


L 


def 




disc 


disc 


L h 


def 


■> H h 


It must be noted however that the previous definition of discrete Hamiltonian system does 
not coincide with the discrete differential embedding of the continuous Hamiltonian system 
which gives 


(A Q) 
(A P) 


dH 



dcy 


(19) 


over T 1 * 1 . We denote by Xa,h the associated discrete vector field 


Xa,h 


dH a dH 
d p°XQ dQ dA ' P ' 


( 20 ) 


We will see that this phenomenon is responsible for the non-coherence between the differential 
and variational discrete embedding of Hamiltonian systems. 


3.3. Discrete Hamilton principle. — A main property of discrete Hamiltonian systems 
is that they are preserving an essential feature of continuous Hamiltonian systems, i.e. the 
variational structure. Precisely, let £a ,h{Q,P) be the discrete functional obtained by the 
A-discrete embedding of the continuous one and defined by 

£ a , h (Q, P) = [Ja(P * A Q- H{Q, P))] n . (21) 

We have : 

Theorem 6. — A couple (Q,P) is a critical point of Ca,h{Q, P) if and only if it satisfies 
the discrete Hamiltonian system associated to H. 

The previous result complete the picture of our construction leading to a full commutative 
diagram between all the object and given by 
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c 


disc 

- > 


A A ,h 


lap 


discrete l.a.p 


Xh 


def 


Xa,v,h 


However, the discrete embedding formalism can be used to obtain a different definition of 
Hamiltonian system using the differential discrete A-embedding. In this case, we obtain a 
non-commutative diagram : 


C 


disc 

- y 


£a ,h 


lap 


discrete l.a.p 


\r disc y / \ ~ 

xH - - > Xa,H T Xa,\7,H 

which is reminiscent of the non-coherence of the procedure between differential and variational 
discrete embedding already in the Lagrangian case (see [6j). 


Proof. The Fredchet derivative of the discrete functional C A h is given for all U, V £ 
Cb(T,K) by 


VC a>h (Q,P)(U,V) = J A (p*AU + V*AQ- dH< ^ P) *U - dH f p P) 

As U £ Co(T, M), we obtain using the discrete integration by parts formula : 
V£ A}H (Q,P)(U,V) = 


J N 


J A (-XP,U + V,AQ- dH ^ P Ku- dH ^ P Kv 


dQ 


dP 


( 22 ) 


J N 

(23) 


which can be rewritten as 


V£ a ,h(Q,P)(U,V) = 


Jv (y * ( AQ _ WQ'P)) _ n Jy P+ mQ,P) 




dP 




dQ 


(24) 


J N 


By definition, a critical point of C a ^h satisfies VC A} h(Q, P)(U, V) = 0 for all U, V G Co(T, M). 
Using the discrete Dubois-Raymond lemma we deduce that 


= m, 

(VP) = -§£, 


over T d 


(25) 

□ 


3.4. Comparison with West and al. and Opris and al. approaches. 
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3. 4- 1- West and al. approach. — In [12j and [8], the authors introduce a definition of discrete 
Hamiltonian systems. The second paper uses the discrete Lagrangian formalism developed in 

m 

The main difference between these approaches and our result is that our definitions are co¬ 
herent with the formalism of discrete embedding developed in [6]. This remark has strong 
consequences. In particular, in the discrete embedding formalism there is no changes between 
the Hamiltonian or Lagrangian defining the continuous system and the discrete analogue. The 
discretisation preserves all the relations between the different objects (Lagrangian, Hamilto¬ 
nian, Legendre transform) globally leading to the previous commutative diagrams. In the 
contrary, in [12] or [8], the discrete Lagrangian associated to a continuous one change the 
Lagrangian function (see |12] .p.363 or [9], p.192-195). Moreover, it breaks the differential 
form of the Euler-Lagrange equation (see §.8.3.5 in [^). 

The same is true in the Hamiltonian case. The authors introduce discrete Hamiltonian func¬ 
tions which do not coincide with the discrete functions associated with the continuous Hamil¬ 
tonian function (see [8],p.5512-5513 and [6], §.3.1). As a consequence, they are able to show 
an analogy between the continuous and the discrete objects (see [12], Part I,§.1.6) but miss 
the fact that this is more than an analogy as proved by the previous diagrams. Using discrete 
embeddings (more precisely the variational one) all the structures and objects are exactly 
transported in the discrete framework. 

3.4- 2. Opris and al. approach. - In (m,§-5), the authors give a definition of discrete Hamil¬ 
tonian systems. This definition is made independently of any discretisation procedure. As 
a consequence, an explicit comparison of their approach is difficult. However, one can find 
some similarities in the definitions even if the global picture given by our framework is dif¬ 
ferent. Indeed, the starting point is a discrete Legendre transform with respect to a discrete 
Lagrangian. As a consequence, they obtain a Hamiltonian function which has the same form 
as in the continuous case (see [7], formula (27) p.26) and a form of the discrete Hamiltonian 
equations which are similar to ours (see [7], formula (32) p.26). However, due to the frame¬ 
work used by the authors, the comparison to the continuous case is very difficult contrary to 
our derivation of discrete Hamiltonian systems. 


PART III 

DISCRETE HELMHOLTZ CONDITIONS 


In this Section, we solve the inverse problem of the discrete calculus of variations in the 
Hamiltonian case. We first recall the usual way to derive the Helmholtz conditions following 
the presentation made by Santilli D3|- We have two main derivations : 

— One is related to the characterization of Hamiltonian systems via the symplectic two- 
differential form and the fact that by duality the associated one-differential form to a 
Hamiltonian vector field is closed. Such conditions are called integrability conditions. 

— The second one use the characterization of Hamiltonian systems via the self adjoint¬ 
ness of the Frechet derivative associated to the differential operator associated to the 
equation. These conditions are usually called Helmholtz conditions. 
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Of course, we have coincidence of the two procedures in the usual case. In the discrete case 
however, it seems that the second one is more appropriate to a generalization, avoiding the 
definition of discrete differential forms as done for example by E.L. Mansfeld and P.E. Hydon 
mm) or Z. Bartosiewicz and al. [3j in the time-scale setting. As a consequence, we 
follow the second way to obtain the discrete analogue of the Helmholtz conditions. 


4. Hemlholtz conditions for Hamiltonian systems 


4.1. Symplectic scalar product. — In this Section we work on M? d ,d > l,d £ N. The 
symplectic scalar product (•, •) j is defined for all 1,1 £ Mi 2d by 


(X,Y)j = (X, JY ), 


where (•, •) denotes the usual scalar product and J = 


0 Id 
-Id 0 


with Id 


on 


We also consider the If 2 symplectic scalar product induced by (■, •) j 
C’HM.R 2 ") by 

= [ (f(t),g(t))jdt . 

J a 


(26) 

the identity matrix 
defined for f,g £ 
(27) 


4.2. Adjoin of a differential operator. — In the following, we consider first order dif¬ 
ferential equations of the form 


fi (q\ = ( X q(QiP)\ 
dt \p) \X p (q,p)J 

The associated differential operator is written as 



q - X q (q,p)\ 
p- X p (q,p)J 


(28) 


(29) 


A natural notion of adjoin for a differential operator is then defined. 


Definition 7. — Let A : C l {\a, b\, M 2ti ) — > C 1 ([a, b], R 2n ) . We define the adjoin Aj 
of A with respect to < • > 1,2 j by 

< A ■ f , g >l 2 ,j=< A*j • 9i f >l 2 ,j ■ (30) 

An operator A will be called self-adjoin if A = Aj with respect to the L 2 symplectic scalar 
product. 


4.3. Hamiltonian Helmholtz conditions. — The Helmholtz’s conditions in the Hamil¬ 
tonian case are given by (see [14], Theorem. 3.12.1, p.176-177) : 

Theorem 8 (Hamiltonian Helmholtz conditions). — Let X(q,p) T = (X q (q,p), X p (q,p)) 
be a vector field. The differential equation T28\) is Hamiltonian if and only the associated 
differential operator given by H29\) has a self adjoin Frechet derivative with respect to the 
symplectic scalar product. 

In this case the Hamiltonian is given by 

H(q,p)= [ [p-X g (\q,\p) - q ■ X p (\q,\p)\d\ 

J 0 


( 31 ) 
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The conditions for the adjointness of the differential operator can be made explicit. They co¬ 
incide with the integrability conditions characterizing the exactness of the one-form associated 
to the vector field by duality (see [14], Thm.2.7.3 p.88). 


Theorem 9 (Integrability conditions). — Let X(q,p) J = (X q (q,p),X p (q,p)) be a vector 
field. The differential operator 0^ b given by [29 1) has a self adjoin Frechet derivative with 
respect to the symplectic scalar product if and only if 

T 


OXg 

dq 


+ 


dX p 

dp 


= 0 , 


dX q 

dp 


dq 


are symmetric. 


(32) 


Of course, the first condition corresponds to the fact that Hamiltonian systems are divergence 
free, i.e. that we have divA = 0. 


5. Discrete Helmholtz’s conditions 

We derive the main result of the paper giving the characterization of first order difference 
systems which ate discrete Hamiltonian systems. 

5.1. Discrete symplectic scalar product. — Following the continuous case, we intro¬ 
duce a discrete analogue of a symplectic scalar product and the notion of self-adjointness for 
discrete difference operators. 


We recall that for X, Y £ Lj the scalar product is given by for all 1,1 6 Lj by 

{X,Y) L 1 A = [J A (X*Y)] N . (33) 

The symplectic version is obtained as in the continuous case. 

Definition 10. - LetX,Y € C'(T,M 2 ' ; ). The symplectic A-scalar product over Lj is defined 

by 

(X,Y) L 2 AJ = (X,J-Y) L , A . (34) 


5.2. Adjoin of a discrete finite difference operator. — We consider first order differ¬ 
ence equations of the form 

aq\ = fx Q (Q,p)\ 

XP) \Xp(Q,P)J- 

The associated finite difference operator is written as 



._(AQ-X q (Q,P) 

I K XP - Xp{Q,P\ 

A notion of symplectic adjoin for discrete finite differences operators can be defined : 


(36) 


Definition 11. Let A: C(T, M 2d ) —» C(T, R 2d ) be a discrete difference operator. 
We define the adjoin Aj of A with respect to (•, ■) L 2 A j by 

{A{X), Y) L 2 A J = (A*j(Y),X) L 2 A}J . (37) 


As in the continuous case, a main role will be played by self-adjoin discrete operators. 






12 


JACKY CRESSON 1 - 2 AND FREDERIC PIERRET 1 


Definition 12. — A discrete operator A is said to be self-adjoin with respect to symplectic 
A-scalar product over Lj if A = A*j. 


5.3. Discrete Helmholtz conditions. — Let X(Q,P) J = (Xq(Q, P),Xp(Q, P)) be a 
discrete vector field. The main result of this Section is the discrete analogue of the Hamiltonian 
Helmholtz conditions for discrete difference equations. 

Theorem 13 (Discrete Hamiltonian Helmholtz conditions) 

Let X(Q,P) J = (Xq(Q, P),Xp(Q, P)) be a discrete vector field. The discrete difference 
equation i35\) is a discrete Hamiltonian equation if and only if the operator 0\ defined by 
Tib]) has a self-adjoin Frechet derivative with respect to the discrete Lj symplectic A-scalar 
product. 

Moreover, in this case the Hamiltonian is given by 

H(Q, P) = [ {P*X q (\Q,\P) -Q*Xp(\Q,\P)]d\. (38) 

Jo 

The proof is given in Section 15.41 


We give an explicit characterization of discrete vector fields satisfying the Hamiltonian 
Helmholtz conditions. By coherence with the continuous case, we call them discrete 
integrability conditions. 


Proposition If (Discrete integrability conditions). The operator 0\ defined by 
f36 1) has a self-adjoin Frechet derivative at ( Q,P ) G C(T, M) x C(T,R) if and only if the 
conditions 


are satisfied over T^. 


dX Q 

dQ 

9Xq 

dP 



0 (CHI) 


and 


dX P 

~0Q 


are symmetric 


(CH2), 


(39) 

(40) 


The proof easily follows from the following Proposition. 


Proposition 15. 

by 


Let U, V £ Co(T,M d ). The Frechet derivative DO(Q, P) of L36}) is given 


DO(Q,P)(U, V) = 


(A U- 


dx Q 


*U- 


ox Q 




dQ *^ dP 


(41) 


and his adjoin DOj(Q,P) with respect to the symplectic A-scalar product is given by 

/ XT /nv- \T 


DOj(Q,P)(U,V) = 


l V W -(W) T * C 7 +(S) T ^ y 


(42) 
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Proof. — Let U, V G Co(T, M d ) and A,Bq C(T, M rf ). The Frechet derivative DO(Q,P) 
follows from simple computations and is given by 


(DO(Q,P)(U,V),(A,B)) L 2 A j = 




J a ( AU*B- i^r*v)*B 


dXc 


-vu*a + 


dQ 


dP 


J N 


Using the discrete integration by part formula, we obtain 


(DO(Q,P)(U,V),(A,B)) L 2 AtJ = 


-u* 

+ u* 


VR- 


A A + 


(dX P \ 

V dQ J 

m 


~k A T 


*A- 


dX, 


Q 


dQ 


dX, 


Q 


dP 


*B 


*B 


N 


As a consequence, the symplectic discrete adjoint of DO is given by 


DO*j(Q, P)(A, B) = 
which concludes the proof. 


(AA + (!&y*A-(!&) 


XB- 


dX P 

dQ 


a A T 


ax Q 

dQ 


*B 




(43) 


□ 


5.4. Proof of Theorem 1131 — If A is a discrete Hamiltonian vector field then there 
exist a function H 6 (7(M d x such that Xq = ^ and Xp = — The discrete 


Hamiltonian Helmholtz conditions of order one are satisfied. Indeed, we have 


d 2 H _ d 2 H 
dPdQ ~ dQdP 


by the Schwarz lemma. 


Reciprocally, we suppose that X satisfies the discrete Hamiltonian Helmholtz conditions. We 
define the function H as 


H(Q, P) = [ [P*A q (AQ,AP) -Q*X P (AQ,AP)]dA (44) 

Jo 

with the functional 

£a,h(Q,P) = [Ja(P * XQ — H(Q, P))] N . (45) 

The Frechet derivative of Ca,h at ( Q,P) along U, V £ Co(T,M) is given by 

VCa,h(Q,P)(U,V) = [J a {P*XU + V*AQ-DH(Q,P)(U,V))\ . (46) 

Using the discrete integration by part formula, we obtain 

VCa.h(Q,P)(U,V) = [Ja(-XP*U + V*AQ-DH(Q,P)(U,V ))] n . (47) 
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We compute DH(Q,P)(U,V) over T^. We denote (*) = ( XQ,XP ). We have 


DH(Q,P)(U,V) = j [V * X Q (*) - U * X P {*) + 
— A Q k 


A Pk 


dx Q ( -). u+ B 2m* v 


9Q 


dP 


BXp{, Ku + ^X*v 


dQ 


dP 


dX 


f dX Q (*) 

V dQ 


+ AU* 


(vkX Q (*) - UkX P (*) + XUk 


kp 


-w-« 


3X q (* 

dP 


dX. 


Using the Hamiltonian Helmholtz conditions we obtain 

r-l 


DH(Q, P)(U, V) = [VkX Q (*)-UkX P {*)-XUk 




dP 


dQ 


+ AU* 




dP 


dQ 


dX. 


Remarking that 


dX P (*) _ dX P (*) dX P {*) l 

<9A <9P + dQ y 

9X Q (*) SX Q (*) „ , 9A Q (*) ^ 

5A “ SP + 5Q ^ 


we deduce that 

/»1 o 

DH(Q, P)(U,V) = j -(XVkX Q (*)-XUkXp(*))dX, 

and finally, integrating with respect to A : 


(48) 

(49) 

(50) 


(51) 


DH{Q, P)(U , V) = V k X q (Q, P) — U k X P (Q , P). (52) 

Replacing DH(Q, P)(U,V) by its expression in the differential of Ca,h, we obtain 

VCaAQ, p )( u , V) = [Ja (V k (A Q - X q (Q, P))-Uk (VP - X P (Q, P)))) N . (53) 

A critical point ( Q,P) satisfies VCa,h(Q, P){U,V) = 0 for all U, V G Co(T, M n ). Using the 
discrete Dubois-Raymond lemma we deduce that 


/ (A Q) =X Q , 
1 (VP) =Ap, 


( 54 ) 


over T^. This concludes the proof. 
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5.5. Comparison with Opris and al. results. — In mm, the authors prove discrete 
analogues of the Helmholtz’s theorem. 

The result in |;I] is closer to our previous work [5] as they consider discrete Lagrangian 
systems. 

In [7], the authors define first discrete Hamiltonian systems and discuss the Helmholtz con¬ 
ditions for first order difference equations (see m,§. e). As already noted in S I3.4.21 it is not 
clear to compare their results with the continuous case. Indeed, the objects are not related 
to discretisation of the continuous one. As a consequence, the problem is in one hand more 
general but in the other one easier to solve. Indeed, no constraints are associated with an 
underlying continuous structure to respect. We can remark that they use also the self-adjoin 
characterization of the Hamiltonian structure (see [7].Prop.9) and not a characterization using 
discrete analogue of one-form and closeness. 


PART IV 
APPLICATIONS 


In this Section, we study several discrete finite differences systems for which we decide it they 
admit a discrete Hamiltonian structure. 


6. The linear case 


Let us consider the following discrete linear system defined over by 


f A Q = aQ + /3P, 
X VP =1 Q + 5P. 


(55) 


where a, /3, 7 and 5 are constants and Q,P € C(T,M n ). The Helmholtz condition (CH2) 
is clearly satisfied. However, the system (15511 satisfies the condition (CHI) if and only if 
a T b — 0. 

As a consequence, Hamiltonian linear difference equations are of the form 


J A Q = aQ + /3P, 
\ VP = 7 q-«p, 


(56) 


Using formula (1381) we compute explicitly the Hamiltonian which is given by 

H(Q, p ) = \ (/3P 2 - lQ 2 ) +&Q ■ P . 


(57) 


7. Newton’s equation 


The Newton’s equation (see [2]) is given by 


j q = p/m , 

1 P = -U'{q ), 


( 58 ) 
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with m G M + and q,p G M. d . This equation possesses a natural Hamiltonian structure with a 
Hamiltonian given by 

H{q,p) = ^—p 2 + U(q). (59) 

2m 

Using the construction of Section 3.1, a natural discretisation is given by 


f A Q = P/m, 

I VP = —U'(Q), 


(60) 


defined over T 1 * 1 . The Hamiltonian Helmholtz conditions are clearly satisfied. 


Remark 2. — It must be noted that the Hamiltonian associated to w is given by 

H(Q,P) = ^-P 2 + U(Q), (61) 

2m 

which can be recovered by formula (1381) . 


As shown in [6] and Section 3.3, there is non-coherence between the differential and variational 
discrete embedding. Moreover we can detail the differences and consequences between these 
two approaches. Consider the A differential embedding of (1581) defined over T + which is given 
by 

A Q = P/m, 

A P = —U'(Q), 

where Q,P G C(T,M). Introducing Z = p(P ) we obtain the following equivalent system over 
given by 

A Q = p(Z)/m, 

VZ =-U\Q). 

Now the problem is that the discrete field is depending on Z through the shift operator p. 
As p(Z) = Z — hV(Z), inserting the definition of p(Z) in (f63l) we obtain 


(62) 


(63) 


/ A Q =1-(Z + hU'(Q)), 
\VZ= -U'(Q). 


(64) 


In that case, the Helmholtz condition (CH2) is clearly satisfied. The condition (CHI) is 
equivalent to U"(Q ) = 0 which means that (1641) is Hamiltonian if and only if U is linear in Q 
which is a very restrictive condition. 


8. Linear friction 

Let us consider the continuous system with friction 

« = J,/m ' (65) 

p = - lp - q , 

where 7 G K and g,p G M. d . As in the previous example, two choices are possible. The one 
mimicking the canonical form of the discrete Hamiltonian systems is given by 

A Q = P/m, 

VP = - 7 P - Q, 

and the one corresponding to discrete differential embedding is given by 

A Q = P/m, 

AP = - 7 P - Q. 


( 66 ) 


( 67 ) 
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In the first case, the Hamiltonian Helmholtz condition (CH2) is clearly satisfied. The condi¬ 
tion (CHI) is equivalent to 7 = 0. In other words, the system (1661) has a discrete Hamiltonian 
structure if and only if there is no friction which is consistent with the continuous one. 


In the second case, using the same trick as in the previous example, we obtain the following 
equivalent system 

A Q = p(Z)/m, 

VZ =- 7 p(Z)-Q, 

where Z = p{P). Assuming that 7 / 1 /h, we obtain 



AQ 

VZ 


= — l-hrjQ- 


(69) 


As previously, there is only one condition to be satisfied which is (CHI). This condition is 
equivalent to 7 = ^. In other words, the system (1661) has a discrete Hamiltonian structure if 
and only if the time step in the discretisation is equal to my. 


9. Modified Harmonic oscillator 


Let us consider the following discrete linear system defined over by 

/ AQ = P + a(P), 

1 VP = Q. 


(70) 


where Q,P G C(T, M 2 ) and a(Q) T = (P2,0). In that case, the Helmholtz condition (CH2) is 
not satisfied. Indeed, is not symmetric because 9 °gp^ 1 = 1 and d °gp^ 2 = 0. Consequently 
the system (1701) does not have a discrete Hamiltonian structure. 


CONCLUSION AND PERSPECTIVES 


As in the Lagrangian case studied in [5], one can generalize our result in the setting of 
the time-scale calculus. Contrary to what happens in the Lagrangian case where technical 
difficulties related to the composition of operators seem to cancel such a generalization, the 
Hamiltonian case is possible due to its linear nature in the discrete differential operators. We 
refer to m for more details. 

Another important task is to defined a notion of discrete differential forms adapted to our 
point of view in order to fully follow the continuous formulation of the integrability conditions 
in term of closed differential forms. 

Last but not least, the previous approach can be used to solve the inverse problem of the cal¬ 
culus of variations in the stochastic setting introduced by J-M. Bismut in [4]. In this book, he 
defined the notion of stochastic Hamiltonian systems and develop the corresponding stochas¬ 
tic calculus of variations providing all the ingredients to generalize Helmholtz’s conditions to 
the stochastic case. 
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